function g=prox_lst(x,a,b,c)

num=size(b,1);
xi=unidrnd(num,1);
dim=length(x);
g=inv(c*a(xi,:)'*a(xi,:)+eye(dim))*(c*b(xi)*a(xi,:)'+x);

